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ABSTRACT 

Context. We present a method for including gas extinction of cosmic-ray-generated UV photons in chemical models of the midplane 
of protoplanetary disks, focusing on its implications on ice formation and chemical evolution. 

Aims. Our goal is to improve on chemical models by treating cosmic rays, the main source of ionization in the midplane of the disk, in a 
way that is consistent with current knowledge of the gas and grain environment present in those regions. We trace the effects of cosmic 
rays by identifying the main chemical reaction channels and also the main contributors to the gas opacity to cosmic-ray-induced UV 
photons. This information is crucial in implementing gas opacities for cosmic-ray-induced reactions in full 2D protoplanetary disk 
models. 

Methods. We considered time-dependent chemical models within the range 1-10 AU in the midplane of a T Tauri disk. The extinction 
of cosmic-ray-induced UV photons by gaseous species was included in the calculation of photorates at each timestep. We integrated 
the ionization and dissociation cross sections of all atoms/molecules over the cosmic-ray-induced UV emission spectrum of Hi. By 
analyzing the relative contribution of each gas phase species over time, we were able to identify the main contributors to the gas 
opacity in the midplane of protoplanetary disks. 

Results. At 1 AU the gas opacity contributes up to 28.2% of the total opacity, including the dust contribution. At 3-5 AU the gas 
contribution is 14.5% of the total opacity, and at 7-8 AU it reaches a value of 12.2%. As expected, at 10-15 AU freeze-out of species 
causes the gas contribution to the total opacity to be very low (6%). The main contributors to the gas opacity are CO, CO2, S, SiO, 
and Ot. OH also contributes to the gas opacity, but only at 10-15 AU. 

Key words. Astrochemistry; Protoplanetary disks; Molecular processes; (ISM:) cosmic rays 
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The midplane of protoplanetary disks has been considered a 
dead zone, because the lack of a source of ionization prevents 
the development of magneto-rotational instabilities, which 
are thought to drive the accretion process. The midplane of a 
disk corresponding to a Class II source ar ound a T Tauri star 
is op a que to stellar and interstellar U V (Ivan Z adelhoff et al] 
200 It IWoitke et all l2009l) and X-ray dGlassgold et alJ l2007t 
" 1201 1) photons, which corresponds to a region 
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located at z/r<0.05 and l<r<10 AU. However, most regions of 
the midplane are far from being dead zones because of cosmic 
rays, which directly ionize the gas and heat the dust grains. An 
important consequence of the interaction of cosmic rays and 
H2 molecules is the emission of a locally generated UV field 
that can ionize/dissociate species. More specifically, regions 
that have a value for Ay of a few are dominated by cosmic-ray 
processing, since they can penetrate colu mn densities of 2 ~ 
150 g cm" 2 dUmebayashi & Naka nolll981l) . 

Hence, the midplane of the disk between 1 and 10 AU can 
be viewed as a cosmic-ray dominated region, where cosmic 
rays are the main source of ionization and therefore the main 
driver of the chemical evolution. Steady-state chemical models 
applied to these regions cannot fully describe the chemical 
evolution of the midplane, as the chemical rel axation timescale 
can be as long as 10 8 yr ( Woitk e et alJl2009|). which is longer 
than t he lifetime of the disk ( Fedele et al J 120101 : lHaisch et alJ 
I20T0H S emenov &Wiebdl201ll) . For this reason it is necessary 



to approach the study of the chemical evolution of the disk 
from a time-dependent model. One of the main catalysts 
of the cosmic-ray driven chemi stry in these regions is OH 
dChaparro Molano & Kampll2012l) . especially for the formation 
ofCOandH 2 0. 

Current chemical models of protoplanetary disks include 
to a large extent the effects of cosmic rays, whether they 
model the chemistry u s ing steady-state dWoitke et alJ 1 2009: 
IWillacv & Woods! l2009t iGorti et al l [20111 iThi et alJ 1201 lh or 
time-dependen t dVisser et alJ I2009L I2OIU IWalsh et alJ l201fi 
ISemenov et al. 2010) solutions. However, the effects of the 
local gas opacity and grain growth in cosmic -ray-induced UV 
processes are largely overlooked, as the parameters for esti- 
mating cosmic -ray-induced photoionization/dissoci ation rates 
are usually taken from mo lecular cloud literature dGredel et al.l 
ll989tlWoodall et al.ll2007h . 



In molecular clou d modeling dCecchi-Pestellini & Aiellol 

[T99l:IShen etal. . 2004) it is customary to ignore those effects, as 
they are not relevant for the cold, gas-poor environment present 
deep inside molecular clouds and the interstellar medium 
(ISM). Another important factor to consider is grain growth in 
protoplanetary disks, which reduces dust UV opacity compared 
to molecular cloudfl This leads to an enhancement in the 
cosmic -ray-induced UV flux (Chaparr o Molano & Kampl|2012l) 



1 See lD'Alessio et al] (200 1 ) for evidence of grain growth from spec- 
tral energy distributions of protoplanetary disks 
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Table 1. Distance from the star, temperature, and density con- 
ditions corresponding to midplane regions in the protoplanetary 
disk structure in Fig.Q] following the PrqD iMo simulation of a 
passive irradiated disk ( Woit ke et al.l l2009). The stellar parame- 
ters used in this simulation are found in Table|2] 



r (AU) 


r(K) 


«<h> (cm 3 ) 


1 


80 


10 14 


3-4 


65 


10 12 


7-8 


40 


10" 


10-15 


20 


10 10 



with respect to previous ISM based estimates dShen et al.l l2004). 
especially in gas-poor regions. This field, which is enhanced by 
a factor 40 at 10 AU, drives gas phase formation pathways for 
saturated molecules that can later freeze on the dust surface. 
Wherever the physical conditions allow for penetration of cos- 
mic rays in disks, our analysis of cosmic -ray driven chemistry 
applies. Since we aim to perform a quantitative analysis, we 
chose the physical conditions from a particular disk model. 
However, our results do not depend on that particular choice. 

This paper is structured as follows. The strategy to imple- 
ment our model in Sect. 2 is followed by a discussion of the 
physical conditions in gas-rich regions of the disk midplane 
in Sect. 3. In Sect. 4 we discuss our treatment of cosmic-ray- 
induced UV photoprocesses including the effects of gas opacity. 
Sect. 5 deals with the particulars of our chemical model, fol- 
lowed by a summary (Sect. 6) and a discussion (Sect. 7) of our 
results. Finally, the main conclusions from our work are summa- 
rized in Sect. 8. 



2. Methodology 

Our goal is to identify the main chemical contributors to the ex- 
tinction of cosmic -ray-induced UV (CRUV) photons in different 
regions of the midplane of a T Tauri disk. The temperature in 
the disk midplane can be high enough to prevent species from 
freezing onto the surface of grains immediately after their gas 
phase formation. This general freeze-out of chemical species 
is found to happen at distances b eyond approximately 10 AU 
(Chaparro Mol ano & Kampll2012l) . Thus, we chose the 1-10 AU 
range of the disk midplane in this work, becuase the temperature 
and density conditions are ideal for studying regions with very 
different gas compositions. For example, the temperature at 7-8 
AU (see Table [TJ coincides with the onset of thermal desorption 
of CO, which will evaporate from the surface of grains while 
leaving the abundances of other frozen species for the most part 
unchanged. 

The physical input conditions for our chemical evo lution 
models are given by a ProDiMo model (IWoitke et al.ll2009t) . and 
are listed in Tables Q] and [2] We implemented a time-dependent 
calculation of the CRUV photorates that includes the extinction 
provided by the dust and also by gas species, which depends 
on their abundance. For this reason, we recomputed the CRUV 
photorate at each timestep of the simulation to account for the 
changes in gas phase abundances. This extinction was then 
integrated over the wavelength range and emission probability 
of CRUV photons. The emission probability was obtained from 




1AU 3-5AU7-SAU 10- AU 
1 10 

r[AU] 



100 



Fig. 1. Density structure model of a T Tauri disk as a function 
of radial distance from the star and the relative height, following 
the ProDiMo simulation dWoifke et al.ll2009l) . The black ovals 
show the regions of the disk according to Table Q] The relevant 
disk parameters are provided in Table|2] 



Lyman and Werner emission lines of H2 that is excited either 
by direct cosmic-ray interactions or by s econdary electrons 
gener ated in cosmic -ray excitation of H2 dPrasad & Tarafdarl 
Il983h . 

Using this scheme, we can trace species that have a strong 
impact on the absorption of CRUV photons for the midplane of 
T Tauri disks. By taking into account the extinction contribution 
of dust grains, we can compare it to the gas extinction, thus ob- 
taining a time-dependent value for the opacities over the CRUV 
wavelength range. We implemented this method in our time- 
dependent chemical rate equation solver chem_compact, which 
we previously used for studying ice formation i n the comet for- 
mation zone dChaparro Molano & Kampll2012h of a passive T 
Tauri disk with a low accretion rate. This code was benchmarked 
against steady- state chemica l abun dances from the ProDiMo 
simulation from lWoilke et al.1 d2009l) . 

2. 1 . Gas-rich regions 

Cosmic ray-induced UV photons can be absorbed by the 
material in the local environment where they are generated. 
Both gas and dust can absorb these photons and become a 
source of local extinction, but this depends on the local density 
and temperature conditions. For instance, at 10 AU most of the 
material is frozen onto the surface of dust grains, which means 
that the extinction of CRUV depends entirely on the local dust 
properties. In the absence of gas extinction, grain aggregation 
in protoplanetary disks can lead to an enhanced CRUV flux 
dChaparro Molano & Kamp] l2012h . In these gas-poor envi- 
ronments CRUV photoprocesses will not be affected by the 
composition of the gas and the chemistry can be described in a 
fairly straightforward fashion. 

By contrast, in regions closer to the central star the environ- 
ment is heated up and most of the chemical species will stay 
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Table 2. Table of modeling parameters. 



Parameter 


Symbol 


Value 


Stellar mass 


M, 


1M G 


Effective temperature 


T e g 


5770 K 


Stellar luminosity 


L, 


1 L<s 


Disk mass 


M D 


0.01M G 


Inner disk radius 


Rin 


0.5 AU 


Outer disk radius 


Rout 


500 AU 


Gas surface density power law index 


6 


1.5 


Dust-to-gas mass ratio 


PdlPs 


0.01 


Minimum dust grain size 


^min 


0.1 


Maximum dust grain size 


^max 


10 /mi 


Mean molecular weight 


A< 


1.35 


Dust grain size power law index 


P 


3.5 


Dust material mass density 


Pgmd 


2.5 g-cm" 3 


Dust grain albedo (UV) 


0) 


0.57 


Dust opacity (UV) 


*uv 


6.8xl0 3 cm 2 -g _1 


Cosmic ray ionization rate (H2) 


&J 


5 x 10~ 17 s _I 


Number of active layers 


^Lay 


2 


Adsorption site area 


A 


0.0/ x 10 cm 


C adsorption energy 


E i 


630 K 


CO adsorption energy 


i-CO 

E b 


960 K 


C0 2 adsorption energy 


r COo 


2000 K 


CH 3 adsorption energy 


E b ' 


920 K 


CH 4 adsorption energy 




hook 


adsorption energy 


F° 
1 


630 K 


2 adsorption energy 


F - 


960 K 


OH adsorption energy 


77OH 
E b 


1000 K 


H2O adsorption energy 




4800 K 


Si adsorption energy 




2100 K 


SiH adsorption energy 


E-SiH 
E b 


2300 K 


SiO adsorption energy 


E-SiO 

E b 


2800 K 


Fe adsorption energy 


£ A 


3300 K 


Mg adsorption energy 


F Mg 


4200 K 



in the gas phase while leaving significantly reduced layers of 
frozen species. Any change in the chemical composition of the 
gas will either enhance the CRUV field or quench it efficiently, 
depending on the CRUV cross section of the dominant species 
in the gas. Thus, if we aim to understand the complex coupling 
effects between chemistry and CRUV photons, we need to 
fully incorporate the contribution of the gas in the local CRUV 
extinction. 

In Fig. Q] we show the regions under study in a plot of 
the density structu re of the disk obtained using ProDiMo 
dWoitke et alJ l2009h . which provides us with a self-consistent 
hydrostatic structure from which we obtain parameters such 
as temperature, density, and intensity of the local UV field 
compared to the ISM Draine field. By studying regions at 
different distances from the star, we can identify the species that 
are locally dominant in their CRUV opacity. 

The specific temperature and density conditions for each spe- 
cific region of the disk midplane are found in Table Q] Near the 
inner rim (at 1 AU from the star) we focus on the role of CRUV 
photoprocesses in the highly efficient OH forming region near 
the inner rim. Moving farther away from the star, the temper- 
ature and the density decrease, which causes more material to 
freeze onto the surface of dust grains. To understand the role of 



freeze-out, we probed the disk at two more regions: at 3-4 AU 
and 7-8 AU. For the 10 AU region we refer th e reader to part I 
of our paper dChaparro Molano & Kampll2012l) . 

3. Cosmic-ray-induced processes 

Cosmic rays penetrate to the disk midplane predominantly from 
the vertical direction because the column density along all other 
directi ons is too high. From th e analysis of cosmic -ray penetra- 
tion in Semenovetal. (2004), it follows that in the midplane 
of the particular generic T Tauri disk chosen in this study (see 
Table|2]for a list of parameters), at radial distances larger than 1 
A U cosmic rays can penetr ate almost unhindered. Recent work 
bv lPadovani & GalLI d201 ll) on the penetration of cosmic rays in 
molecular clouds suggests that interaction with magnetic fields 
is more important than previously thought. However, those 
results do not necessarily apply for the particular magnetic field 
geometry of protoplanetary disks, and more detailed modeling 
is necessary to clear the picture of cosmic -ray and magnetic 
field interactions. 

The main chemical byproducts of direct cosmic-ray ioniza- 
tion are H+ (from H+) and He + dHerbst & Klempererl fl973h . 
While He + is very good at dissociating molecules and passing 
on its charge to the products of the reaction, H3 hydrogenates 
and ioni zes CH compounds, and helps cre ate water from atomic 
oxygen (Chap arro Molano & Kampll2012l) . 

Cosmic rays can also ionize the medium in a more subtle 
way: by inducing a UV field that comes from secondary ioniza- 
tion of molecular hy drogen. The process, know n as the Prasad 
Tarafdar mechanism (Prasad & Taraf dial [19831) . starts when an 
electron with a typical energy of 30 eV is released after the 
cosmic -ray ionization of an H2 molecule. This secondary elec- 
tron can also ionize another H2 molecule. A value for the total 
rate £h 2 of both direct and secondary ionization of H2 is not en- 
tirely agreed upon, but a conservative value of 



( Hj =5x 10~ 17 s" 



(1) 



has been obtained both for H^ measurements in the ISM 
dlndriolo et al.ll2007h and fro m theoretical estimations based on 
measured cosmic-ra y spectra (Cecchi-Pest ellini & Aiellol ll992: 
iMicelotta et 10120 Ilk . 

The emitted electron then hits another (neutr al) H2 molecule, 
which leaves it in an excited electronic state dSternberg et al.1 
Il987h iRiahi et al. 2006), after which it spontaneously decays to 
the excited vibrational states of the B 'E* and C 1 S l4 levels. In 
the subsequent decay to excited vibrational states of the ground 
electronic level X 'E* of H2, Lyman and Werner photons are 
emitted in the 90-170 nm range. These cosmic-ray-induced UV 
photons can then either ionize/dissociate a gas species or hit a 
dust grain. 

We define the CRUV photoprocess efficiency as the fraction 
of CRUV photons that dissociate a species and are not locally 
absorbed by the gas or dust: 



-I 



3.28 PHz 



P(v)<77(v) 

1.76 PHz °totW 



dv . 



(2) 



Here P is the emission probability profile of a CRUV photon, 
cr, is the photoprocess cross section (in cm 2 per species), and 
<r tot is the total (gas+dust) cross section, which is a measure of 
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the local extinction. If a given species with a high CRUV cross 
section is very abundant in the gas phase, the photo rate of that 
species will have a maximum value and will be low for all other 
species in the gas. This shielding effect cannot be ignored in 
regions where high-density/temperature combinations create a 
gas-rich environment. 

Finally, cosmic rays can directly heat dust grains and cause 
desorption of ices. We took this effect into account, but it is more 
predominant in outer regions of the disk. 

3.1. CRUV emission probability profile 

The emission probability profile of CRUV is obtained from the 
transition probability of the first three electronic levels (B 'Z^ 
and C l T, u ) of molecular hydrogen. The cross section for exci- 
tation of H2 into a level v' is proportional to the optical band 
oscillator strength /o v <. 



oyo K /ov' ■ 



(3) 



We use proportionality here, as it is enough to obtain a normal- 
ized emission probability profile. The oscillator strength is 



/ov ocAv-o — 



(4) 



Here g is the statistical weight. The emission probability for a 
transition from the level v'J' of an electronically excited state ; 
to the v'J' level of the ground electronic state is then proportional 
to the Einstein A coefficient for emission and the cross section 
for excitation: 



(5) 



This probability is normalized over all transitions to v" J" levels 
of the ground electronic state: 



f' A ' 

JOV v'J'V'J" 



Tjv'J' Tjv"J" /fjV^W.v" 



(6) 



Each transition probability is then convolved into a Voigt line 
profile to account for both natural and thermal broadening: 



/ , (v;v^)= J P(v^)0 v (v-v;, J ,) 



(7) 



The emission probability profile for each excited electronic state 
is then 

" (8) 



Thus, the probability that a CRUV photon will be emitted after 
an H2 cosmic -ray ionization process is 



P(v) = P B '^( V ) + P c '^(v) + P c '^(v) (9) 



This probability is normalized over the frequency range 1 .76 - 
3.28 PHz (90-170 nm). We obtained the Einstein A coefficients 
and the frequencies for the relevant transitions from the tables of 
Abgr all et alj (2000), and the statistica l weights were calculate d 
from the guidelines in Appendix A of Och kin & Kittelll ([2009). 



3.2. CRUV gas opacity 

The total cross section cr tot in Eq. © is (in cm 2 per hydrogen 
atom) 

cr tot (v) = 5-^(1 - co) + J] &<7v(v) . (10) 

j 

Here oj is the grain albedo, = h ; /h<h> is the abundance of 
the species /', and cr^' is the grain UV extinction cross sec- 
tion per hydrogen atom. We ref er the reader to Eq. (10) in 
Chaparro Mol ano & Kampl d2012l) to see how the cross section 
relates to the dust UV opacity and other dust parameters in Table 
[2] Since the dust UV extinction curve is fairly flat in the CRUV 
frequency range, we used a frequency average for this value. Due 
to photon conservation, the previous expression leads to the fol- 
lowing relation: 



^(1-^)7^ + 6^ = 1. 



Here y dust is 



Tdust 



P(y) 

CTtot(v) 



dv . 



(11) 



(12) 



Thus measures the fractional contribution from the species 
i to the CRUV extinction. The shape of the radiation field that 
dissociates or ionizes a species F(v) not only depends on the 
CRUV emission probability, but also on the CRUV extinction of 
all other species cr tot (v): 



F(y) 



P(v) 
0" to t(v) 



(13) 



The frequency-dependent opacity for a species i can be written 

as 

n/u\ 

(14) 



Ki(v) = & o-j(v) — cm 2 g (gas) . 

Pgas 



Here p g as=«(H>/U'«H, with n being the mean molecular weight 
(see Table|2]i. We can now define a "gray" (frequency averaged) 
opacity as 



<*/> = 



r 3.28 PHz 

Ji.76phz F(yMv)<*y 

r 3.28 PHz 

Jl.76PHz F ^ dV 



(15) 



Using Eqs. ( fT3T > and ([T4l this expression can be rewritten as 



<*;> 



&7i 



TdustjU'WH 



(16) 



With this expression, we can accurately measure the opacity of 
each species throughout a time-dependent chemistry run. The 
variables y, and ydust are calculated at every timestep of the 
simulation as part of the CRUV photoprocess rate calculation. 
Given that the CRUV photoprocess rate constant k, oc y u the re- 
sulting rate is then coupled to the abundance of all other species 
that contribute to the opacity. 

Since photoionization/dissociation of species can be contin- 
uum and/or line processes, we convolved the line photoprocess 
cross sections with a Voigt profile in order for both natural and 
Doppler line broadening. The total cross section is then the sum 
of the line and continuum processes for each molecule. Values 
for the cross sections were obtained from the tables in the 
Leiden photoprocess database hosted by E. van Dishoeck 
http : //www . strw . leid enuniv . nl/ " ewine/photo/ 
( Ivan D ishoec k et al.ll2006h . 
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Table 3. Table of chemical species in the chemical networks. # 
indicates an ice species. 



Type 


Symbol 


Atoms 


H, He, C, O, S, Si, Mg, Fe 


Ions 


He + , Si + , Fe + , H~, H + , C + , 




+ , S + , Mg + 


Molecules 


H 2 , H 2 0, CH 2 , HCO, SiO, C0 2 , 




SiH, CH 3 , CH 4 , OH, 2 , CO, CH, H 2 CO 


Molecular 


HCO + , CH+, H+, SiH + , SiO + , 


Ions 


CH+, H 3 + , H 3 + , SiH J, CH+, 




CH+, H 2 + , SiOH + , CH + , H+, 




OJ, CO + , OH + , CO+ 


Ice 


C#, CO#, C0 2 #, CH 3 #, CH 4 #, 




0#, 2 #, OH#, H 2 0#, 




Si#, SiH#, SiO#, Fe#, Mg# 



It is safe to consider only the effects of dust CRUV ex- 
tinction if we know a priori that most of the gas species are 
frozen onto the surface of grains. This means that if gas phase 
abundances are low enough to be neglected in Eq. ( fTOb , dust 
grains will be the sole contributor to opacity. This approach is 
frequently taken in models that include CRUV photoprocesses 
(Ster nberg et all 1 19871: Idredel etal.lfl987l 119891: IWoodall et all 
2007), mostly because it is valid for ISM and molecular cloud 
conditions (low molecular gas abundances). 

At 10 AU the environment is cold enough (r gas =20 K, 
"<h>=10 10 cm -3 ) to study the chemical evolution without takin g 
into account the gas opacity dChaparro Molano & Kamp]|2012f) . 
Under these conditions, the CRUV rate constant takes a simpli- 
fied form that does not depend on species abundances. 

4. Chemical model 

The code chem_compact, de s cribed in our previous pa- 
per dChaparro Molano & Kampl 1201 2ft . is our VODE based 
(Bro wn et alj [l989) gas/grain chemical rate equation solver. In 
it we include a rea ction network based on the Umist06 database 
for astroc hemistry (W oodall et al. 2007) including H2 formation 
on grains (ICazaux & Tielensl l2002) and ad/desorption reactions: 
Adsorption a n d thermal and stell a r UV p hotode sorpti o n from 
lAikawa et all (fT996h : iLeeeretalJ (fl985ft: Obe rg et all (12009). 
cosmic-ray direct desorption from lHasegawa & Herbsj (11 993) 
and c osmic-ray-induced photodesorption from iRoberts et alJ 
Surface reactions are not considered because they are 
beyond the scope of this paper. 

Table[3]lists all gas and ice species considered in our model. 
In this work w e use the low metal initial abundances from 
IJenkinsI (|2009h . where absorption lines of these metals are 
measurecQ from various clouds against a bright background star. 
The low abundances imply that Si, Fe, and Mg condense into 
dust grains before the formation of the disk. This is particularly 
relevant for the formation of SiO and its maximum abundance 
levels, although formation of CO and H2O is not affected. 



2 Despite assuming significant metal deplet ion, the metal abundan ces 
here are about a factor 10 higher than in Graedel et al] dl982l) or 
iLee et 



Table 4. Table of initial abundances, atomic (column 2) and 
molecular (column 3). # indicates an ice species. 



Species 


Abundance 


Species 


Abundance 


(Atomic) 




(Molecular) 




H 


1 


H 2 


5.00x10-' 


He 


8.51xl(T 2 


He 


8.51xl0~ 2 


O 


3.31xl(T 4 


CO 


1.76xl0 4 


C 


1.78xl(T 4 


H 2 0# 


1.09xl0 4 


Fc 


5.01x10-" 


O 


3.55xl0" 5 


Mg 


2.24x10-" 


OH 


5.66x10" 


S 


1.90x10-" 


Fe# 


5.01x10-" 


Si 


1.78xlQ-" 


H 


2.56xl0~ 6 






H,0 


2.55x10-" 






Mg# 


2.24x10"" 






S 


1.73x10-" 






H 2 CO# 


1.63x10-" 






SiO# 


1.02x10-" 



Initial conditions for our disk model were obtained by 
running our chemical evolution code under molecular cloud 
conditions (T = 20 K, «(h> = 10 6 cm 4 ) from atomic low- 
metallicity abundances. The resulting abundances after 10 7 
yr were used as initial conditions. Table [4] lists the initial 
abundances for the molecular cloud run (atomic) and those for 
the disk model. Atomic abundances (column 2 in Table |4]i were 
used as input for the molecular cloud run, which yielded the 
molecular abundances (column 4 in Table [4]i that we took as 
initial abundances for our disk model. 

We decided to ignore sulphur chemistry beyond ionization 
of atomic S because its chemical network is only weakly 
coupled to CO f ormation and is decou pled from other species 
considered here (Sternberg & Dalgarnol ll995ft . Atomic sulphur 
was included as it contributes to the metallicity, and its ioniza- 
tion can be an important source of charge exchange for other 
species. By ignoring other S-bearing species, we can still arrive 
at an upper limit for the importance of atomic S in the CRUV 
gas opacity. 

To estimate the effects of gas opacity in the CRUV photopro- 
cess rates and therefore in the ch emical evolution, we replaced 
the relevant rates from the umist (IWoodall et aLll2007l) database 
with our own rates (Sect. |3). As these rates depend on we re- 
calculated them at every timestep, using the new values for the 
gas abundances. This means that at every timestep, the integral 
in Eq. (O needs to be evaluated from the previous timestep. 

5. Results 

To estimate the relevance of gas opacity for protoplanetary 
disks, we used an estimate of the gas ionization a nd photodis- 
sociation cross sections from iGredel et alJ (1 19891) to compute 
the mean photo dissociation cross section for various species 
in the wavelength range of CRUV photons (90-170 nm). Using 
this value, we computed the abundance that is required for an 
individual species to produce a gas opacity equal to the dust 
opacity. We did this for both molecular cloud and protoplanetary 
disk dust properties. 

H2 does not contribute to the opacity because most of it is in 
the fundamental vibrational level of the ground electronic state 
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Table 5. Comparison between the gas and dust opacity:average 
photoprocess cross section (column 2), maximum gas phase 
abundance required to produce an opacity equal to the dust opac- 
ity for ISM type dust (column 3) and for protoplanetary disk type 
dust (column 4). 



Species cr, (cm 2 ) 


Req. «,/w <H > 


Req. n,/n (H > 


(i) (90-170 nm) 


(PPD) 


(ISM) 





J.OJA 1U 


18 


T. OS V 1 f)~ 
J.70A 1U 


-5 


S ?f)y 1 fr 


-4 


V 1 1 


U.ULA1U 


8 


2.31x10" 


-5 


1 fSRylfV 

1.UOA1U 


-5 


v 1 1 


4 S.fSxlfr 


8 


J. 1JA1U 


-5 


4.12x10- 


-4 


TH, 
v 11 


1.JUA1U 


7 




-5 


9.47x10" 


-5 


V I 1 : 


1 .JLA 1U 


6 


1 . 1 UA 1U 


-6 


1 ^?yl(T 

1.JZA1U 


-5 




U.1JA1U 


6 


L.JUA 1U 


-7 


J.iUA 1U 


-6 


OH 


7.82x10" 


X 


1.96x10" 


-5 


2.56x10- 


-4 


H 2 


4.00x10" 


X 


3.83x10- 


-5 


2.87x10- 


-4 


Fe 


1.80x10" 


s 


8.49x10" 


-5 


1.11x10 


-3 


Mg 


9.20x10" 


9 


1.66x10' 


-4 


4.12x10- 


-4 


CO 


1.82x10" 


7 


8.40x10" 


-6 


1.10x10- 


-4 


Si 


2.14x10" 


7 


7.14x10- 


-6 


9.33x10- 


-5 


HCO 


2.99x10" 


6 


5.12x10- 


-7 


6.69x10 


-6 


SiH 


6.35x10" 


7 


2.41x10- 


-6 


3.15x10- 


-5 


H 2 CO 


3.58x10" 


6 


4.27x10- 


-7 


5.59x10- 


-6 


o 2 


7.20x10" 


X 


2.13x10- 


-5 


5.77x10- 


-5 


S 


3.36x10" 


X 


4.55x10- 


-5 


5.95x10- 


-4 


co 2 


9.10x10" 


x 


1.68x10" 


-5 


5.81x10- 


-5 


SiO 


2.07x10" 


6 


7.41x10- 


-7 


9.68x10 


-6 



level through two main pathways. It can be photodissociated by 
CRUV into C and O, which later recombine in the following 

way: 

CO- > CO; C -A CH 2 -> CO . (17) 

CO is mainly formed via electron recombination of HCO + . 
However, SiO also acts as a catalyst for the formation of CO 
because its reaction with HCO + is very rapid. Hence, whenever 
SiO is present, it efficiently converts HCO into CO, without SiO 
having to be very abundant itself. The OH to SiO to CO chemical 
pathway is the following: 



OH 2 



HCO+ 



SiO CO,SiOH + ; SiOH + Si,OH/SiO . 

(18) 

HCO + is formed in the reaction of CO with H^, which as noted 
above in Section [3] is a byproduct of cosmic-ray ionization of 
H 2 . This cycle ensures a steady supply of OH as a primer for 
sustained formation of SiO and CO, as seen in Fig.|2^b. 
CH4 is also steadily formed from He + ionization of CO: 



CO^C+^ 



+ CO 
CI~^2— >3— >5 * CH4 



(19) 



although at a very low abundance (10 7 ). The latter part of that 
reaction is curbed because the first step is reversed by 2 : 



C + + 2 — > CO + o + . 



(20) 



Atomic oxygen is depleted after 10 4 yr, favoring the formation of 
2 and C0 2 . Abundances of C0 2 along with OH are enhanced 
due to the reactions 



X Therefore it is unable to absorb CRUV that are generated 
in transitions from excited electronic states to vibrationally 
excited levels of the X 'Z^ state. 

Using the results from Table|5] we can check to which extent 
the simplifying assumption holds that gas opacity is negligible 
compared to that of the dust for CRUV photoprocesses. At re- 
gions located at distances of r<8 AU chemical abundances of 
the listed molecules easily exceed the boundary values of Table 
[5] This means that the effects of gas opacity have to be included 
in chemical models of the midplane of protoplanetary disks in- 
side 10 AU. 



CO 



Ht 



HCO + 



Fe 



HCO -> OH,CO/C0 2 



(21) 



The OH enhancement is also curbed by the very efficient 2 /C0 2 
formation reactions, which help keep the atomix oxygen abun- 
dance low: 

OH + O -> 2 + H , (22) 



OH + CO -> C0 2 + H 



(23) 



Atomic oxygen is constantly generated thanks to the steady 
CRUV photodissociation of CO, 2 , and C0 2 . This atomic oxy- 
gen partly goes into the formation of water vapor, which is 
promptly adsorbed onto the grain surface due to the temperature 
and density conditions at 1 AU: 



5.1. Chemistry at 1 AU 

In this region, the moderately high temperature (80 K) en- 
sures that most of the chemical species stay in the gas phase. 
Adsorption of gas phase H 2 is much more efficient than 
thermal desorption, which causes gas phase H 2 to have a very 
low abundance. CO, 2 , C0 2 , SiO, CH4, and oxygen in atomic 
form are the dominant species in the gas phase, whereas H 2 
and to a lesser extent SiO are the only species with a significant 
ice abundance. Fig. [2^ shows the time-dependent evolution 
of the species abundances when considering an appropriate 
protoplanetary disk dust size distributiorQ and including the 
effects of CRUV gas opacity. 

CO is very abundant in the initial conditions run, and re- 
mains the most abundant species besides H 2 . CO is kept at this 



3 Our choice of dust size distribution parameters is typical of a Class 
II source, but limited to current knowledge based on observations. 



O -A OH + ,H 2 + ; OH + — * H 2 + — * H 3 + — » H 2 0/OH . 

(24) 

This water vapor formation pathway explains the long-term abil- 
ity of the system to keep forming water ice, through which it be- 
comes the main oxygen reservoir at very long timescales. OH is 
recovered partly from SiOH + and also from the following HsO + 
reaction, which can form water or OH at the same branching 
ratio: 

H 3 + H 2 0, SiOH + . (25) 

Combining reactions (fT8l and d25l > we can see how near the end, 
the oxygen in 2 transfers to H 2 : 



o 2 



H,0+ 

SiO — -> H 2 . 



(26) 



Gas-phase water created in the previous reaction freezes out al- 
most immediately, as mentioned above (see Fig. |2^). In Fig. [3] 
we can see that after 10 4 yr the gas opacities reach a maximum 
of about 40% with respect to the dust opacity, which remains 
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1 1 1 r 

(a) PPD dust with gas opacity 



(b) ISM dust, no CRUV gas opacity 



CO 

co 2 

CH A 



2 
SiO 
O 



H 2 0# 
SiO# 
CH3 



(c) PPD dust, no CRUV gas opacity 



_l_ 



3.5 



4.5 5 5.5 
log(time/yr) 



6.5 



Fig. 2. Chemical abundance evolution at a distance of 1 AU from 
the central star: (a) Using an ISM dust value for <x^ and ignor- 
ing CRUV gas opacity, (b) using an appropriate <x^ for proto- 



planetary disk conditions and ignoring CRUV gas opacity, and 
(c) using an appropriate <x^ for protoplanetary disk conditions 
and including the effects of CRUV gas opacity. 




3 4 
log(time/yr) 



CO 



SiO 



CO, 



Fig. 3. CRUV gas opacity vs. time at 1 AU. For reference, the 
dust opacity is 68 cm 2 /g (gas). 



Table 6. Ratios between various species abundances and CO at 
1 AU after 10 6 yr for three cases: Including the effects of gas 
opacity, considering only dust extinction from ISM dust param- 
eters, and considering only dust extinction from protoplanetary 
disk dust parameters. 



Ratio 


With gas 


No gas opacity 




opacity, cr™ yr 


ISM PPD 

U dust " dust 


H 2 0#/CO 


0.69 


1 0.66 


o 2 /co 


0.09 


0.19 0.28 


C0 2 /CO 


0.01 


0.33 0.07 


CH 4 /CO 


7x10"* 


2x1 0~ 3 10~ 5 



5.1 .1 . Comparison to our previous work 

To understand the effects of gas opacity, we also ran the simu- 
lation considering dust as the only sink for CRUV photons. We 
used the simplified formula for CR UV photoprocess rates from 
IChaparro Molano & Kampl d2012l) . where the total extinction 
cross section is simply the grain extinction cross section. Since 
this calculation is not dependent on the gas phase abundances, 
the rate coefficients are constant throughout the simulation. 



constanfl. The late O2 enhancement (Fig. |2^) implies that it is 
providing most of the CRUV gas extinction, and consequently 
SiO absorbs fewer CRUV photons. Thus O2 slows down atomic 
oxygen formation via photodissociation of SiO, which is a more 
efficient formation mechanism than the CRUV photodissocia- 
tion of O2. This causes a rebound effect after 10 s yr, as O2 needs 
a constant supply of atomic oxygen to form. Because OH is free 
to react with CO (instead of reacting with O), CO2 formation 
is enhanced. However, CRUV photodissociated atomic oxygen 
from SiO and CO is necessary to sustain CO2 production via re- 
action G.W . This causes CO2 to reach a saturation level near the 
end of the simulation. 



4 For simplicity, we did not include self-consistent grain growth mod- 
els here. 



When considering a UV extinction cross section 
based on a protoplanetary disk dust size distribution, the 
CRUV photon flux is higher th an for ISM dust parameters 
( IChaparro Molano & Kampl |2012|) . This is because dust coag- 
ulation processes in protoplanetary disks increase the average 
size of dust grains, which in turn decreases the effective area 
for dust absorption of UV photons. For protoplanetary disk 
conditions the UV extinction cross section crj^ is 13.5 times 

lower than the typical ISM value, cr^ = 2 x 1(T 21 crrT 1 . 

Fig.[2p shows the time-dependent evolution of the chemical 
abundances for the ISM-like dust parameters used adopted for 
the calculation of CRUV rates in Umist. Here, the main carriers 
of carbon are CO and CO2, followed by the slightly less abun- 
dant CH4 and CH3. After a few 10 5 years the carbon is almost 
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Fig. 4. Chemical abundance evolution at a distance of 3-5 AU 
from the central star. 
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Fig. 5. CRUV gas opacity vs. time at a distance of 3-5 AU. For 
reference, the dust opacity is 68 cm 2 /g (gas). 



equally divided among CH4, CO, and CO2. SiO is the second 
most important gas phase oxygen carrier up to a few Myr, when 
it stops to form. This is because the grains are cool enough to 
hold a large amount of water ice on their surface, effectively 
stopping OH formation and the SiO formation pathway ( fT8l ). As 
the oxygen in CO, CO2, and SiO is transfered to water ice near 
the end of the disk lifetime, carbon binds with HJ and H2, which 
enhances the formation rates of molecules such as CH3 and CH4. 

When comparing the chemical evolution obtained when 
using ISM dust parameters (Fig. \2fa) vs. appropriate protoplan- 
etary disk dust parameters (Fig. [2};), the impact of the larger 
CRUV flux on the abundance of various species is evidenced. 
As mentioned above, CRUV photodissociation processes are up 
to 13.5 times more efficient in protoplanetary disks than when 
estimated using the interstellar UV extinction cros s section 
value used in Umist (Chaparr o Molano & Kampl l2012'). Because 
the effects of gas extinction are neglected, the CRUV flux is 
enhanced uniformly for all species. This affects specifically the 
long-term formation of SiO, O2, and hydrocarbons. For instance, 
CH3 and CH4 are not efficiently formed because they are very 
sensitive to CRUV photoprocesses. CH4 photodissociates into 
CH 2 , and CH 3 photodissociates into CH, CH 2 and CH+ . If the 
CRUV field is strong enough, CH3 and CH4 will be destroyed 
at the expense of CO on a timescale similar to their formation 
timescale. 

The CRUV photodissociation of CO, SiO, and CO2 is very 
efficient for high CRUV fluxes. These photoprocesses steadily 
produce atomic oxygen at a higher rate than when considering 
the effects of gas opacity. This efficient formation of atomic 
oxygen forms OH from HCO in reaction ( l2TT i and preserves the 
OH, SiO cycle in reaction 11 1 St . This means that SiO is never 
depleted even at extremely long timescales (below 10 8 yr). 
Even though CO2 is also created from OH at the expense of 
CO, it never reaches a high abundance because it is destroyed 
via CRUV photodissociation, which ensures that the CO2 
abundances stay low compared to those seen in Fig.|2^. 

Most notably, the CO2/CO, O2/CO, and CO/CH 3 , 4 ratios 
change enormously between the three cases, as shown in Table 



[6] This shows that the chemical evolution obtained when consid- 
ering the effects of gas opacity cannot be obtained by a simple 
interpolation between a low and a high CRUV field, but it has to 
be studied for each species separately. 

5.2. Chemistry at 3-5 AU 

Figure [4] shows the chemical evolution in a region located at 3- 
5 AU from the cent ral star, where «<h>=10 12 cirT 3 and T=65 K 
(Woi tke et al.l2009t) . In this region, reactions (Qj]l ,<0, d26l and 
(l24l still take place. However, as the temperature here is lower 
than at 1 AU, SiO is now mostly in the ice phase, and retains a 
small gas phase abundance. It is kept at a stable but low level 
throughout the simulation by the following feedback reaction, 
similar to reaction ( TT8b : 

SiO ^ SiO + ^ SiOH + SiO/Si, OH . (27) 

Therfore, OH formation is not as efficient as it was at 1 AU. 
The low OH formation rate cannot sustain the formation of O2 
in reaction d22l) when it is forming H2O, and the oxygen in O2 
transfers to H2O via reaction [26] As the O2 abundance starts to 
decrease, the backwards reaction ( l20l i is stopped, which allows 
CO to form CH4 more efficiently via reaction [19] (see Fig. |4j. 
The consequent decrease in O2 opacity seen in Fig.[5]causes the 
CRUV photorates to increase for other molecules, as it allows 
more CRUV photons to dissociate or ionize other species. This 
is shown in Fig [5] where the SiO opacity increases as the O2 
opacity decreases. Now CH4 CRUV photodissociation becomes 
very efficient at forming CH2, which causes the following reac- 
tion 

CH 2 + O -> CO + H + H , (28) 

to be favored over reaction ( l22l . This exacerbates the depletion 
of O2, as seen at the end of the simulation in Fig. [4] The gas 
opacity in this region (Fig. O is less than half the value at 1 AU 
(Fig. [3] note the different vertical scale), which is caused by the 
lack of efficient O2 formation and the freeze-out of SiO. The first 
part of reaction (l2"TT i is reversed by the following reaction: 

HCO + + e - -> CO + H , (29) 
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Fig. 6. Chemical abundance evolution at a distance of 7-8 AU Fig. 7. CRUV gas opacity vs. time at a distance of 7-8 AU. For 



from the central star. 



reference, the dust opacity is 68 cm /g (gas). 



which keeps the CO2 levels low. This is enabled by an increase 
in effectivity of electron recombination reactions vs. neutral re- 
actions brought about by the low-density conditions, compared 
to the 1 AU region. 

5.3. Chemistry at 7-8 AU 

In this re gion, the physical conditions are «<h>=10 h cm -3 and 
T-50 K (Woit ke et alJ l2009). As we can see from the chemical 
evolution plot in Fig. [6] H2O and SiO are completely frozen af- 
ter only 1 yr, with negligible gas phase abundances. However, 
the chemical evolution for the other significant species remains 
more or less the same. The lack of SiO in the gas phase impacts 
the late formation of H2O that appears at 1-5 AU. Since O2 is 
not consumed via reaction d26l its depletion is not as dramatic 
as in the 3-5 AU region. Now the main reaction that removes O2 
from the gas is 

C + 2 -> CO + O . (30) 

This reaction, which does not have an activ ation energy barrie r 
and is very rapid at very low temperatures (ISmith et al.l 12004). 
allows reaction (T% to form CH4 very late in the simulation in 
the same way that we observed for the previous region (3-5 AU). 
Reactions (I17I22I24I I are also responsible for the formation of 
CO, O2, and H2O, respectively. Just as in the 3-5 AU region, 
CO2 is not efficiently formed because its formation via reaction 
(|2TI) is stopped by the backreaction ( 129b . The gas opacity here 
(Fig- ED is ver y similar to that observed at 3-5 AU (Pig. [5j. The 
absence of SiO opacity due to freeze-out enhances the CRUV 
rates for all other species, including CO. 

6. Discussion 

We have included the gas opacity of all atoms and molecules 
that are photoprocessed by cosmic -ray-induced UV photons in 
our chemical evolution model. By doing this, we were able to 
assess whether a given species is relevant for the extinction of 
CRUV in the midplane of a protoplanetary disks. By studying 
the impact that the gas opacity of individual species has on the 
chemistry and vice versa, we can simplify the chemical model 
by including only the most important species. As some opacities 



increase, an important shielding effect appears and changes the 
CRUV photoprocess rates. This affects the chemical formation 
pathways for both gas- and ice-phase species. 

As discussed in Section |4] the most important species 
that contribute to the gas opacity are O2, SiO, S, CO, and 
CO2 consistently throughout the 1-10 AU region of the disk 
midplane. We found that the opacities of O2, SiO (and CO2 in 
the warm regions around 1 AU) are the most variable in time 
and hence are the ones that require most attention. As the CO 
abundance remains more or less constant in those regions, its 
opacity remains fairly constant as well. Even though we did 
not include sulphur chemistry beyond ionization, its CRUV 
opacity ranges from 2.5 to 3.5 cm 2 /g gas (see Figs. [3] EJ [7] and 
[8), reaching its highest contribution with respect to the total gas 
opacity (70%) and with respect to the total dust and gas opacity 
(3.4%) at 10-15 AU. 

When comparing the gas opacity to the total CRUV opacity 
in Table|7l it is evident that its effect cannot be easily dismissed. 
In the region 1-2 AU the contribution from the gas is up to 30% 
of the combined gas and dust opacity. At 7-8 AU, the contribu- 
tion is up to 10%, mostly because of SiO freeze-out. Finally, at 
10 AU the temperature becomes sufficiently low to freeze CO, 
CO2 and CH4 onto the grain surface, and the contribution of the 
gas opacity becomes very low (5%). Therefore, future chemical 
models that aim to include this treatment of CRUV gas and dust 
opacity can consider the contribution of the most significant 
species that our models yield for the probed regions of the disk. 

The general chemical composition of the disk midplane ob- 
tained using our models compare favorably to that in Wals h et all 
(12010 "). particularly in the freeze-out of water at distances above 
1 AU, and in the significant gas phase presence of CO (and CO2 
to a lesser degree) throughout the studied region. We notice the 
same low abundances of H2CO and OH. These similarities are 
caused by the temperature structure, which controls the adsorp- 
tion and desorption processes. 
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Table 7. Total gas opacity as a function of time and distance. 
Percentage relative to the total dust (68 cm 2 /g gas) and gas opac- 
ity is given in parentheses. 



Distance 


Total gas opacity <K 6as ) (cm 2 /g gas) 




10 4 yr 


10 5 yr 


10 6 yr 


1 AU 


23.6 (25.8%) 


25.4 (27.2%) 


26.7 (28.2%) 


3-5 AU 


11.5 (14.5%) 


11.1 (14%) 


9.2(12%) 


7-8 AU 


9.5 (12.3%) 


9(11.7%) 


5.4 (7.4%) 


10-15 AU 


3.6 (5%) 


3.7 (5.2%) 


4.2 (5.8%) 



6.1. The role of OH 

The availability of OH caused by the high and steady abundance 
of SiO at 1 AU causes a late O and O2 enhancement. The 
products of SiO dissociation later form H30 + , which via 
dissociative recombinations (Si, e-) forms water vapor. Atomic 
oxygen is very abundant at 1 AU thanks to cosmic -ray-induced 
dissociation of CO and SiO. However, oxygen-driven OH 
formation acts as a catalyst for steady CO formation. A rebound 
effect is observed at a late stage (10 s yr), as O2 self-shields 
the CRUV photons that drive the OH driven pathways for its 
formation. 

Even though SiO is formed in reaction (fT8l at 3-5 AU, un- 
der these conditions it is very efficiently adsorbed onto the grain 
surface. This affects the formation of OH and thus O2 is de- 
pleted at long timescales, favoring CO formation via CH4 and 
later CH2 reactions with O. As the gas phase SiO is depleted, no 
O2 is transformed into H2O, which avoids the complete deple- 
tion of O2. At 7-8 AU there is an overall enhancement of CRUV 
rates because of this. After 10 AU, the total gas opacity is suffi- 
ciently low for OH to drive the chemistry. Indeed, Fig. [8] shows 
the gas opacity at 10 AU, in which OH is seen to contribute. 
This agrees with the results from our previous research paper 
(IChaparro Molano & Kampll2012l) . 

6.2. Survival of SiO in theAy- 1 region 

The longevity of SiO in our models of the disk midplane (see 
Figs. [2k and |4]i could have potential consequences beyond the 
midplane of the disk. We explored the chemistry at 1 AU in a 
direction perpendicular to the disk plane upwards to densities of 
«<h> = 10 11 ctrT 3 . We saw that the abundance of SiO remains 
steadily high, meaning that SiO can be formed even above the 
midplane. 



Following Aikawal (120071) . the vertical drift timescale r v d of 
a molecule formed in the midplane of the disk that moves up to a 
distance z = A above the midplane can be estimated via the disk 
viscosit>0 v (Shakura & Sunv aevlll973h : 



A 1 

Tvd = • 

V 



(31) 



The value of v = ac s h depends on the sound speed c s , the scale 
height h of the disk and the Shakura Sunyaev a parameter, for which 
we chose a value of 0.1 corresponding to a young disk with M si 
10~ 8 M s /yr. This mass accretion rate v alue suggests an age fo r the disk 
in the range 10 5 5 -10 65 yr, according to Hart mann et al.1 ( I1998h . 




Fig. 8. CRUV gas opacity vs. time at a distance of 10 AU. For 
reference, the dust opacity is 68 cm 2 /g gas). 



A simple calculation for A = 0.1 AU, corresponding to the Ay 
1 region (at r — 1 AU) of the disk, yields 



r vd ^ 280 yr . 



(32) 



If the SiO photodissociation timescale at Ay — 1 is longer than 
the vertical drift timescale, it should be possible for SiO to ac- 
cumulate on the surface where it could be observed. We can 
roughly estimate this timescale, which depends on the strength 
of the lo cal UV field y. Usi ng the detailed radiative transfer 
model of IWoitke et al. (2009), we derive a value for x m th e 
range 1-10 at the Ay — 1 line. Accordingly, an SiO molecule 
can survive for approximately 10 3 yr if \ — F This means that 
it may be possible for SiO to accumulate around Ay — 1 and 
be detectable. For example, SiO masers are detected in the cir- 
cumstellar regions around asymptotic giant branch (AGB) stars, 
and their maser lines are found to originate f rom regions be- 
low the dust evaporation radius (lHabinglll996l) . More recently, 
gas ph ase SiO has been detected in the debris disk around r\ 
Corvi dLisse et al.l 120121) . Hence it should be possible to de- 
tect accumulated SiO even if there are dust grains present. Our 
low ISM metal abundances already take Si incorporation into 
dust grains into account. Given that the region where SiO is be- 
ing formed has a small radial extent (r < 1-5 AU), detec- 
tion of SiO could prove to be very difficult in the near future, 
except for high-sensitivity, high-spatial-resolution observations 
that ALMA could carry out. 



7. Conclusions 

We have developed a tool for including the effects of gas 
extinction of cosmic -ray-induced UV photons in the chemical 
evolution in the disk midplane. The most important species 
that contribute to the gas opacity in the 1-8 AU region of the 
disk midplane are O2, SiO, S, CO, and CO2. Because the CO 
abundance is steady in those regions, its opacity remains fairly 
constant as well. Even though CO competes with water ice for 
being the most abundant species, CO is not a major contributor 
to the opacity, whereas CO2 play a more important role in the 
opacity while not being as chemically abundant as CO. 
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We find that the opacities of O2, SiO (and CO2 in the warm 
regions around 1 AU) are highly variable in time, and their 
effects on the gas and ice chemistry cannot be easily dismissed. 
Therefore, future models should include the gas opacity of these 
species in the calculation of CRUV photoprocess rates. 

As most of the e mission in recent IR observations of OH 
and H 2 Q in disks dCarr & Najital l2008t iNaiita et alj 1201 Ot 
iFedele et alJ 120111: iMeeus et alj 120121) is likely coming from 
layers well above the midplane, the catalytic role of OH in 
the midplane is not necessarily something we can prove by 
observing its surface abundance. 

Finally, we found that if SiO is steadily created in the mid- 
plane at long timescales, it can drift vertically upwards and ac- 
cumulate around the Ay — 1 region, which suggests a possibility 
for future detections of SiO in protoplanetary disks. 

Appendix A: Cross section data 



Table A.l. Cosmic ray-induced UV photoprocess list. 



Species 


Proc. 


Products 


Source 


Type 


C 


Ion. 


C++e- 


Leiden 


Cont. 


CH 


Dis. 


C+H 


Leiden 


Line/Cont. 


CH + 


Dis. 


C + +H 


Leiden 


Line/Cont. 


CH 2 


Ion. 


CH++e- 


UMIST06 


Cont. 








(Gredel et al. 1989) 




CH 2 


Dis. 


CH+H 


Leiden 


Line/Cont. 


CH 3 


Ion. 


CH++e- 


UMIST06 


Cont. 








(No ref . data) 




CH 3 


Dis. 


CH 2 +H 


Leiden 


Line 


CH 3 


Dis. 


CH+Ht 


Leiden 


Line 


CH 4 


Dis. 


CH 2 +H 2 


Leiden 


Line/Cont. 


OH 


Dis. 


O+H 


Leiden 


Line/Cont. 


H 2 


Dis. 


OH+H 


Leiden 


Cont. 


Mg 


Ion. 


Mg + +e- 


Leiden 


Line/Cont. 


CO 


Dis. 


C+O 


Leiden 


Line 


Si 


Ion. 


Si + +e- 


Leiden 


Cont. 


HCO 


Dis. 


CO+H 


Leiden 


Line 


HCO 


Ion. 


HCO + +e- 


UMIST06 


Cont. 








(Gredel et al. 1989) 




SiH 


Dis. 


Si+H 


Leiden 


Line 


H 2 CO 


Dis. 


CO+H, 


Leiden 


Line/Cont. 


o 2 


Dis. 


O+O 


Leiden 


Line/Cont. 


o 2 


Ion. 


OJ+e- 


Leiden 


Line 


S 


Ion. 


S + +e- 


Leiden 


Cont. 


co 2 


Dis. 


co+o 


Leiden 


Line/Cont. 


SiO 


Dis. 


Si+O 


Leiden 


Line 


Fe 


Ion. 


Fe + +e- 


Leiden 


Cont. 



Table I A. 1 1 presents the information on all CRUV photo- 
processes cojisid^rjedherejbased on the UMIST06 reaction 
database (IWoodall et al.1120071) . Column 2 lists whether the pho- 
toprocess ionizes or dissociates the species. Column 3 shows 
the products of the reaction. Column 4 references the source 
for the photoprocess cross section, which for most of them 
is Ewine van Dishoeck's Leiden photoprocess database 
http : //www . strw . leid enuniv . nl/ ewine/photo/ 
( Ivan Dishoeck et ail 120061) . Column 5 refers to whether the 
photoprocess is a discrete and/or continuum absorption process. 

Acknowledgements. We would like to thank Malcolm Walmsley for asking the 
basic question that triggered this research. 
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